Electronic Transport in Dual-gated Bilayer Graphene at Large Displacement Fields 
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We study the electronic transport properties of dual-gated bilayer graphene devices. We focus 
on the regime of low temperatures and high electric displacement fields, where we observe a clear 
exponential dependence of the resistance as a function of displacement field and density, accom- 
panied by a strong non-linear behavior in the transport characteristics. The effective transport 
gap is typically two orders of magnitude smaller than the optical band gaps reported by infrared 
spectroscopy studies. Detailed temperature dependence measurements shed light on the different 
transport mechanisms in different temperature regimes. 

PACS numbers: 72.80.Vp, 72.20.-i, 73.20.Hb, 73.22.Pr 



The ability to electrostatically tune and deplete the 
charge density in two-dimensional electron gases enables 
the fabrication of basic mesoscopic devices, such as quan- 
tum point contacts or quantum dots, which enhance 
our understanding of electronic transport in nanostruc- 
tures [1]. Creating such electrically tunable nanostruc- 
tures in monolayer graphene, a novel two-dimensional 
system is far more challenging due to its gapless na- 
ture. In this respect, Bernal-stacked bilayer graphene 
(BLG) is an interesting material, because of the possi- 
bility of opening a band gap by breaking the symmetry 
between the top and bottom graphene sheets [3|-[5j . 

The low-energy band structure of free-standing BLG 
is gapless but in the presence of an on-site energy dif- 
ference between the bottom and top layers a band gap 
develops. Different methods have been employed to in- 
duce a band gap including molecular doping, coupling to 
the substrate, and electric displacement field generated 
by gate electrodes [6l-ll2|. However, the low-temperature 
(< 100 K) transport characteristics of dual-gated BLG 
devices do not exhibit the strong suppression of conduc- 
tance expected given the large band gaps (up to 250 meV) 
measured by infrared spectroscopy [9l4l2j. In addition, 
only very weak non-linearities were found in the current 
versus source-drain voltage (/-Vsd) characteristics [9j, 
in contrast with the strong non-linear behavior of typ- 
ical semiconducting devices. A more complete study is 
needed to address the transport characteristics of gapped 
BLG devices as well as the role played by disorder. 

In this letter, we report on the electronic trans- 
port properties of dual-gated (back-gated (BG) and top- 
gated (TG)) BLG devices. We focus on the regime of 
large transverse electric displacement fields, 0.8V/nm 
< \D\ < 2.5 V/nm, over 3 times larger than in previous 
low-temperature experiments [9]. Upon the application 
of a large displacement field, we observe an exponential 
dependence of the device resistance on \D\ and density, 
which is accompanied by a strong non-linear behavior in 
the I-Vsb characteristics. However, the size of the ef- 
fective transport gap is on the order of a few meV, two 
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FIG. 1: (a) Optical image of a BLG (outlined by white line). 
A, B, C, and D are contact electrodes and TG is a 1 /im- 
wide top gate electrode. The red and black diagrams are 
set-ups for two and four probe measurements respectively, 
(b) Schematic diagram of the measured device (not drawn 
to scale), (c) When BLG is subject to a transverse electric 
field, a potential difference is induced between top and bottom 
layers, (d) Band structure of free standing BLG (dashed lines) 
and band structure of BLG subject to transverse electric field 
(solid lines). 



orders of magnitude smaller than the optical band gaps 



at the same D [ll|, ll2| , suggesting a strong role played 
by disorder. Temperature dependent measurements in 
the 300 mK to 100 K range show that the conductivity 
follows an activated behavior with three different activa- 
tion energies, including a nearest neighbor hoping regime 
at the lower temperatures. However, it is the conduction 
mechanism at intermediate temperatures (2-70 K) which 
is responsible for most of the temperature variation of 
the conductivity of our devices, in which conductivity 
increases by several orders of magnitude. 

BLG devices [Fig. []Ja)] are fabricated by mechanical 
exfoliation [13|, followed by the deposition of Cr/Au con- 
tact electrodes, AI2O3 growth, and deposition of a top 



gate [14[ . The typical mobility of our devices before oxide 
growth is between 1500 and 2000 cm 2 /V-s. The mobility 
can degrade significantly after oxide deposition (by ~30% 
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FIG. 2: (a) Differential resistance at zero bias as a function 
of top and back gate voltages at T = 300 mK (log scale). 
The black horizontal (dashed) and black diagonal (solid) lines 
correspond to zero charge densities in non top-gated and top- 
gated regions, respectively, (b) Cuts in (a) at different dis- 
placement fields D (colored lines) and at n = (black line). 
Each cut corresponds to the lines in (a) with the same color, 
(c) 3D plot of (a), (d) Schematic band structure and Ey. 
As D is varied at fixed n, the size of the band gap changes 
while Ey remains fixed. We can also shift Ey and keep the 
band gap constant by varying n and keeping D the same, (e) 
I-Vsb characteristics at different values of D for n = at 
T — 300 mK. (f) The resistance as a function of carrier den- 
sity and offset top gate voltages at various D's. The black 
dashed lines are linear fits in log-scale, (g) The slope a of the 
linear fits in (f) as a function of D. The slope a decreases 
linearly with decreasing D. 



in the device shown) which indicates that additional im- 
purities have been introduced to the system during the 
oxide deposition. Charge and resonant impurities have 
been proposed to be dominant sources of scattering in 
graphene [15l-ll9j . The impurities alter the potential pro- 
file experienced by the charge carriers and lead to the 
formation of electron and hole puddles I?], Eo|, El | . They 
also induce tail states as well as localized states inside the 
band gap [22|-[25j. In addition, in BLG devices, charged 
impurities could lead to a spatial variation of the band 
gap. This suggests that adsorbates and/or charged impu- 
rities in the oxide may play a significant role in electronic 
transport in BLG. 

We first focus on the transport properties at 300 mK. 
Fig. [2(a) shows the zero-bias resistance measured be- 



tween electrodes A and B in Fig. QJa). A parallel plate 
capacitor model yields a charge density under the top- 
gated region n = C B g(Vbg - Vbg) + Ctg(V T g ~ V£ G ) 
where C is the capacitive coupling, V is the gate volt- 
age, and (Vq G , Vj? g ) is the charge neutrality point (CNP) 
in the top-gated region. Following the convention from 
Zhang et al [11[, we define the average electric displace- 
ment field D — (.Dbg + D*tg)/2. I n this letter, we con- 
centrate on the case when .Dbg = D*tg = £ bg(%g — 
^bg/^bg) where 6bg = 3.9 is the relative dielectric con- 
stant of Si02 and ^bg = 285 nm. The sharp rise in re- 
sistance along the n = [Fig. [2(b)] is characteristic of 
BLG (9I. [Tol. [26j .The maximum on-off ratio we can achieve 
in this device at 300 mK is on the order of 10 6 , with a 
minimum resistance of ~300 Q [measured in a four probe 
geometry at (Vbg, Vtg) = (-170, -10) V]. We note that 
our high quality oxides enable us to apply a displacement 
field over 3 times larger than in previous low temperature 
experiments [9|, resulting in an on-off ratio and insulat- 
ing resistivity well over three orders of magnitude larger 
at 300 mK. Such insulating behavior makes BLG a good 
candidate for the fabrication of electrostatically designed 
mesoscopic devices. 

The black diagonal curve in Fig. [2(a-b) shows a slice 
of the resistance for n = 0. Along this curve, we vary 
D in the range [— 2.5, +1.8] V/nm, while keeping Ep in 
the top-gated region at the CNP [top path in Fig. 12(d)]. 
Beyond \D\ ~ 1 V/nm the resistance exhibits a clear ex- 
ponential behavior with increasing \D\. Such resistance, 
if arising from a thermally activated behavior across a 
band gap, would be proportional to exp(E g (D) /2&bT) 
where E g (D) is the .D-dependent band gap, ks is Boltz- 
mann's constant, and T is the temperature of the system. 
However, a fit to our data using the E g (D) obtained by 
infrared spectroscopy [11], yields an effective T w 70 K, 
which is much higher than the 300 mK at which the mea- 
surement is performed. Hence, the conduction we observe 
cannot be explained by activated behavior across such an 
optical band gap: disorder plays an important role, and 
the associated energy scale for transport is about two 
orders of magnitude smaller. This exponential increase 
in resistance is also accompanied by the development of 
strong non-linear transport characteristics. Figure [2(e) 
shows measurements of the DC current between A and 
B electrodes as a function of bias voltage and D, while 
keeping n = 0. We observe a clear non-linear behavior, 
which is consistent with the presence of an effective trans- 
port gap. However, the onset of non-linearity occurs on a 
scale of a few meV (2.2 meV for \D\ ~ 2.5 V/nm), again 
about two orders of magnitude smaller than the optical 
band gap. 

We now consider the behavior of the resistance ver- 
sus gate voltage at various constant D [Fig. [2(b), colored 
lines]. By sweeping the gates along a constant D line, we 
effectively hold the size of the band gap fixed and shift the 
Fermi energy Ey from the valence band to the conduction 
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FIG. 3: (a) Conductivity as a function of inverse temper- 
ature from 4.2 - 100 K. (b) Conductivity as a function of 
inverse temperature from 300 mK to 100 K. The black and 
red curves are the fits to the equation [1] with NNH and VRH 
terms respectively, (c)-(f) The extracted parameters from the 
fits plotted as a function of \D\. 



band [bottom path in Fig.EJd)]. We observe an exponen- 
tial decrease in resistance as we sweep the gate voltages 
away from the CNP (n = 0), followed by a slower de- 
crease, which we associate with Ep reaching the valence 
and conduction band mobility edges. The decrease ap- 
pears symmetric on both the electron and hole sides, and 
depends on D. Figure [2ff) shows a few traces in detail, 
together with the fitted straight lines in the exponential 
regions, and Fig. [2jg) shows the slope a of these lines as 
a function of D, which exhibitis an approximately linear 
behavior. 

To gain further insight into the transport properties, 
we study the temperature dependence of the conductiv- 
ity in the 300 mK to 100 K range at different D's, and 
at n = 0. Figure [3fa) shows the conductivity measured 
using a four probe geometry in the temperature range 
4-100 K for nine different values of D, while Fig. [3jb) 
shows the conductivity over the entire temperature range 
for five different values of D. In contrast to previous low 
temperature experiments at low D-fields [9[, we observe 
two distinct temperature regimes, which become more 
pronounced at high D: a two-component fast decrease 
of the conductivity from 100 to 2 K, followed by a rela- 
tively weak T-dependence from 2 K to 300 mK, which we 
associate with thermally activated hopping (see below). 

At intermediate to high temperatures, we find that 



the conductivity is well described by an activated be- 
havior with two different activation energies, 1/p = 
Y^i=i l/Pi Gxp(Ei/kBT) where p and E are resistiv- 
ity and activation energy respectively [black lines in 
Fig. [3fa)]. The extracted fit parameters (pi,p2,^i, and 
E2) are plotted in Fig. [3jc-d). The higher energy scale, 
Ei, may be related to thermal activation across the op- 
tical band gap [lOj. However, due to the very limited 
high-T range (70-100K), our data barely probe the onset 
of this exponential increase, and the values of E\ ob- 
tained are very likely underestimated. The intermedi- 
ate energy scale, E2, exhibits a linear dependence on D, 
and is therefore approximately proportional to the band 
gap [ll|. However, its value is two orders of magnitude 
smaller than the observed optical band gap at the same 
\D\ (HI- Still, it is the conduction mechanism in this in- 
termediate T-regime that is responsible for most of the 
measured variation in conductivity of the devices. 

The overall behavior of the conductivity with T is 
reminiscent of that observed in disordered semiconduc- 
tors [27], where transport via impurity bands and ther- 
mally activated hopping dominate transport at low tem- 
peratures. However, BLG devices are unique, in that 
their band gap can be continuously tuned by electrostatic 
means, and also because its two-dimensionality means 
that all disorder effects are surface rather than bulk ef- 
fects. Compressibility measurements 28], Hij have shown 



that a very large density of states exists in gapped BLG 
on Si02, even in moderate D fields. Most of these states, 
however, are localized and do not contribute to trans- 
port, and in general the relationship between density of 
states and conductivity is more complex in the insulating 
regime. 

To explain the origin and some of the qualitative fea- 
tures of the conductivity at intermediate T, we consider 
a model of an impurity band arising from the interac- 
tion of negatively charged donors [30]. Such model has 
been used to explain the observed conduction in the in- 
termediate temperature regime in germanium semicon- 
ductors [HI, [HI . Disorder can lead to the formation of 
an impurity band in which charge carriers are localized 
and a hopping mechanism dominates the conduction as 
we will show below. At low temperature, the major- 
ity of these localized states are empty or singly occupied. 
However, some of the states can become doubly occupied 
which leads to an extra band at higher energy due to the 
Coulomb interaction. Carriers in this band are weakly 
localized and hence possess a higher mobility. We can es- 
timate the localization length of these states by equating 
E2 with a charging energy e 2 /er where e ~ 5 is the aver- 
age dielectric constant of our top and bottom oxides, and 
r is the localization length. This estimate yields a local- 
ization length on the order of 100 nm (at \D\ = 2.5 V/nm) 
to 1 /im at (\D\ = 0.9 V/nm), comparable with the width 
of the top gate. Such large localization length supports 
that the states are weakly bound. While this model ac- 
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counts for some of the features present in the data, a 
rigorous theoretical model which takes into account the 
geometry and particularities of disorder in BLG, beyond 
the scope of this paper, is needed for a direct compar- 
ison and understanding of the complex relationship be- 
tween transport and density of states, specially in the 
non-linear regime. In addition, compressibility measure- 
ments at larger D and lower temperatures [H, Q may 
have enough resolution to observe the effects of disorder 
at these energy scales. 

The conductivity between 2 K and 300 mK decreases 
weakly with temperature, which indicates that we enter 
a hopping conduction regime through strongly localized 
states [27j. We perform a fit to the conductivity for the 
complete temperature range and all D with 



P piexp(^) P2 exp(^_) S' 

where H = psex.p(Es/kBT) for nearest neighbor hop- 
ping (NNH) or p'z exp(Tp/T) 1//3 for variable range hop- 
ping (VRH) [27]. Both the NNH and VRH fits agree 
reasonably well with the data [Fig. [3jb)]. VRH has 
been proposed to be the transport mechanism for gapped 
BLG at intermediate and low T [9[. However our mea- 
surements yield a value of p f 3 that grows exponentially 
with \D\ [Fig. [3fe) inset]. Such strong dependence is un- 
expected because, in VRH, the factor exp^o/T) 1 / 3 al- 
ready includes the strong exponential contributions from 
both hopping between sites and differences in energy lev- 
els [I?}. Hence, we propose that the transport mechanism 
in our BLG devices is NNH in the temperature and D 
regime explored. 

For NNH, ps = p® exp(2r / a) where r is the dis- 
tance between hopping sites and a is the localization 
length [27(1 . The magnitude of r can be approximated 

— 1/2 

from the density of impurities, n^, by r = n i . The 
CNP in this device is located at (Vg G ,VP G ) « (16,0) V, 
which corresponds to m ~ 10 12 cm -2 and r « 10 nm. 
The linear fit from Fig. [3^e) yields a localization length 
a ~ 4/[D(V/nm)]nm. This allows us to estimate the 
crossover temperature T VRH at which the conduction 
mechanism changes from NNH at high T to VRH at 
low T. As the temperature is lowered, it is feasible 
for electrons to hop to further sites but closer in energy 
due to the reduced coulomb interaction. This transition 
takes place when 2r/a and E^/fe^T are comparable [27j, 
which yields T VRH < 80 mK for \D\ < 2.5V/nm. This 
crossover temperature is almost 4 times smaller than the 
lowest T we have studied (300 mK) which supports that 
NNH dominates the conduction. 

In addition, the hopping activation energy E3 de- 
creases with decreasing \D\ [Fig. Eff)]. This is an indi- 
cation of the BLG making a transition from a strong to 
a weakly insulating state. As the band gap gets smaller, 
the electron-hole puddles start to merge and hence cre- 
ate a channel for carriers to percolate from one electrode 



to the other [33]. Extrapolating the fit to £3 = 0, we 
obtain a value for the displacement field corresponding 
to this transition of \D\ = 0.7V/nm corresponding to 
Eg/2 « 35 meV which is consistent with BLG Ep fluctu- 
ations for an impurity density rii ~ 10 12 cm -2 . 
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